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Abstract 

Let P be a polygonal curve in M d of length n, and S be 
a point-set of size k. We consider the problem of find- 
ing a polygonal curve Q on S such that all points in S 
are visited and the Frechet distance from P is less than 
a given e. We show that this problem is NP-complete, 
regardless of whether or not points from S are allowed 
be visited more than once. However, we also show that 
if the problem instance satisfies certain restrictions, the 
problem is polynomial-time solvable, and we briefly out- 
line an algorithm that computes Q. 

1 Introduction 

Measuring the similarity between two geometric objects 
is a fundamental problem in many fields of science and 
engineering. However, to perform such comparisons, a 
good metric is required to formalize the intuitive con- 
cept of "similarity." Among the many metrics that have 
been considered, Frechet distance has emerged as a pop- 
ular and powerful choice, especially when the geometric 
objects are curves. Shape matching with Frechet dis- 
tance has been applied in many different fields, includ- 
ing handwriting recognition [7' , protein structure align- 
ment [5], and vehicle tracking [3]. 

In this paper, we consider the basic problem of mea- 
suring the similarity of two polygonal curves. However, 
in our problem, the input is only partially defined. In- 
stead of being given both curves, we are given only one 
polygonal curve P as well as a point set S. Our problem 
is to complete this partial input by constructing a polyg- 
onal curve Q that best matches the given curve, under 
the restriction that the constructed curve's vertices are 
exactly S. 

Since our metric of choice is Frechet distance, we be- 
gin with an popular, intuitive description the concept. 
The metaphor of a person walking a dog is often used, 
with the dog walking along one curve and its owner walk- 
ing along the other. The Frechet distance between the 
two curves is the length of the smallest leash that would 
allow both the person and the dog to reach the end of 
their respective curves without ever backtracking or let- 
ting go of the leash. If a very short leash is sufficient, 
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Figure 1: A problem instance and its solution. 



then the curves are very similar. But if a longer leash 
is required, then the curves are very different. 

Suppose our dog owner wants to walk his dog while 
walking down a path in a given park. The extremely 
curious and territorial dog wants to sniff and/or mark 
every single tree along the path, running directly from 
tree to tree. Mindful of their dog's proclivities, the 
owner goes shopping for a leash long enough to allow 
the dog to have its way without pulling the owner off 
the path. For a given path and set of trees, will a leash 
of a given length be sufficient? More formally, given a 
polygonal curve P, a point-set S, and a real number 
£ > 0, does there exist a polygonal curve on S that vis- 
its every point in S and has Frechet distance less than 
e from PI Unfortunately for owners of territorial dogs, 
we show in this paper this problem is NP-complcte. 

2 Previous Work and New Results 

The decision version of the Frechet distance problem 
asks, given two geometric objects and a real number e > 
0, is the Frechet distance 6f between the two objects 
less than el Alt and Godau p] showed that, when the 
objects in question are polygonal curves of length n and 
m, this problem can be solved in 0(nm) time. They also 
showed that finding the exact Frechet distance between 
the two curves can be done in 0(nmlog(nm)) time. 

Maheshwari et al. [5] examined the following variant 
of the Frechet distance problem, which we refer to as the 
Curve/Point Set Matching (CPSM) problem. Given a 
polygonal curve P of length n, a point set S of size k, 
and a number e > 0, determine whether there exists a 
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polygonal curve Q on a subset of the points of S such 
that Sf(P, Q) < £• They gave an algorithm that decides 
this problem in time 0(nk 2 ). They also showed that the 
curve of minimal Frechet distance can be computed in 
time 0(nk 2 \og(nk)) using parametric search. 

Wylie and Zhu [5] also explored the CPSM problem 
from the perspective of discrete Frechet distance. In 
contrast to the continuous Frechet distance, which takes 
into account the distance at all points along both curves, 
the discrete Frechet distance only takes into account the 
distance at the vertices along the curves. They formu- 
lated four versions of the CPSM problem depending on 
whether or not points in S were allowed to be visited 
more than once (Unique vs. Non-unique) and whether 
or not Q was required to visit all points in S at least 
once (All-Points vs. Subset) They showed that, under 
the discrete Frechet distance metric, both non-unique 
versions were solvable in 0(nk) time, and both unique 
versions were NP-complete. 

In this paper, we show that the Continuous All-Points 
versions of the CPSM problem, both Unique and Non- 
unique, are NP-complete. Table [1] shows the eight ver- 
sions of the problem, with our results highlighted. 





Discrete 
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NP-C [8] 


Open 
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All-Pts Unique 


NP-C [5] 


NP-C* 


Non-Unique 


P [8] 


NP-C* 



Table 1: Eight versions of the CPSM problem and their 
complexity classes. New results starred. 



3 Preliminaries 

Given two curves P 7 Q : [0,1] M d , the Frechet 
distance between P and Q is defined as Sf(P,Q) = 
inf ff , T max te[0 ,i] ||P(cr(t)), Q(r(t))\\, where cr,r : [0, 1] -> 
[0, 1] range over all continuous non-decreasing surjective 
functions [4]. We make use of two commonly noted ob- 
servations: 

Observation 1 Given four points a^b, c, d G M. d , if 
\\ac\\ < e and \\bd\\ < e, then Jf(o0, cd) < e. 

Observation 2 Let P\, P2, Qi, and Q2 be four curves 
in R d with 5 F (Pi,Qi) < £ and 5 F {P 2l Q2) < e. If the 
ending point of Pi (resp. Q\) is the same as the starting 
point of P2 (resp. Q2) then 6f(Pi + P%, Qi + Q2) < £, 
where + denotes concatenation. 

We now give a number of geometric definitions, some 
of which were used in [BJ. For a given a point p £ M. d 
and a real number e > 0, let B(p, e) = {q € M. d : \\pq\\ < 



e denote the ball of radius e centered at p, where ||-|| 
denotes Euclidean distance. For a line segment L C M. d , 
let C(L, s) = Upei ^(p, e) denote the cylinder of radius 
e around L. Note that a necessary condition for two 
polygonal curves P and Q to have Frechet distance less 
than £ is that the vertices of Q must all lie within the 
cylinder of some segment of P. 

Let the continuous function P : [0, 1] — > R d represent 
a curve in M. d . Given two points u,v G P, we use the 
notation u -< v if u occurs before v on a traversal of P. 
The relation >- is defined analogously. 

4 Restricted Satisfiability Problem 

Our NP-completeness result is obtained via reduction 
from a restricted version of the well-known 3SAT prob- 
lem. The 3SAT problem takes as input a boolean for- 
mula with clauses of size 3, and asks whether there exists 
an assignment to the variables that makes the formula 
evaluate to TRUE. If we restrict the input to formulas 
in which each literal occurs exactly twice, the problem 
becomes the (3,B2)-SAT problem. This may seem to 
be a rather extreme restriction, and, indeed, formulas 
of this type with less than 20 clauses are always satisfi- 
able. However, despite this restriction, the problem was 
shown to be NP-complete in [2], and an example of an 
unsatisfiable formula with 20 clauses was presented. 

In order to simplify our reduction, we make the fur- 
ther restriction that no two clauses have two literals in 
common. In other words, we restrict the input to formu- 
las in which the function that maps each literal to the 
pair of clauses it appears in is injective. For any formula 
that violates this assumption, an equivalent, compliant 
formula can easily be constructed using the "balanced 
enforcers" described in [2]. We therefore assume formu- 
las to have this property for the remainder of the paper. 

In the following sections, we first give a brief summary 
of the construction. We then describe the the main 
gadget used, building it incrementally Finally, we give 
the full construction in detail. 

5 The Reduction 

Let $ be a formula given as input to the (3,B2)-SAT 
problem. We construct a polygonal curve P and a point 
set S such that $ is satisfiable if and only if there exists 
polygonal curve Q whose vertices are exactly S with 
Frechet distance less than e from P. 

First, we construct a gadget consisting of components 
of P and S that will force any algorithm to choose be- 
tween two possible polygonal path constructions. The 
gadget is constructed in such a way that these two 
choices are the only possible polygonal paths along the 
gadget's component of S with Frechet distance less than 
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(a) The two dashed curves are the only (b) More corners have been added. (c) All except one of the two interior 

possible curves on 5* with Frechet dis- The two interior points are only reach- points are covered, regardless of which 

tance less than 1 from the given (solid) able on different curve possibilities. path is chosen, 
curve. 



Figure 2: The separation gadget, step by step. 



e from P. These two path possibilities will correspond 
to TRUE and FALSE assignments for a given variable. 

Then, we create a series of points in S to represent 
the clauses in $, one point for each clause. For each 
variable, a gadget will be placed so that the pair clause 
points representing the clauses in which the variable's 
positive instances occur are only reachable along one of 
the two curve possibilities, and likewise for the negative 
instances. Once this has been done for each variable in 
any polygonal curve Q who vertices are exactly S 
with Sf(P,Q) < £ will correspond to an assignment to 
the variables of $ in which every clause is satisfied, thus 
making the formula evaluate to TRUE. Furthermore, if 
no such curve exists, then there can be no such satisfying 
assignment for 

5.1 Separation Gadget 

We begin the description of our main gadget with a spe- 
cific example, which we will proceed to generalize. Con- 
sider the problem instance shown in Figure I2al with 
S = {a,b,c,d}. It is clear that the answer to this 
problem instance is "no" ; no polygonal curve on S with 
Sf(P,Q) < £ can visit both b and c. However, suppose 
this P and S were part of a larger problem instance. 
Suppose further that other segments of P come within 
£ of b and c. The answer to the problem instance is 
no longer so obvious. Even if both points cannot be 
reached the first time they are encountered, it is possi- 
ble that whichever point was skipped could be covered 
in the future. This creates the fundamental difficulty 
that leads to our reduction. 

Figure I2bl shows an extension of the previous config- 
uration, with more corners, all symmetrically the same 
as the first. Note how we have not increased the num- 
ber of options; there are still only two possible paths to 



take. We can add as many of these corners as we like 
without breaking this property, as long as they all bend 
in the same direction. 

The corner points must be placed very precisely to en- 
sure the above properties hold. Because their position is 
so constrained, using them to represent elements of $ in 
our construction would be difficult. At each corner, the 
two path possibilities alternate between the boundary 
of the cylinders and the interior. As shown in Figure 
I2b[ extra points in the cylinder interior are still only 
visible from the other interior points, and therefore we 
can add as many as we like without affecting the path 
possibilities. Thus, by extending the segments between 
the corners, we can create large regions which are only 
reachable along one of the two possibilities. 

There is still a problem to be addressed; as more and 
more corners are added, more and more points are cre- 
ated that would be skipped by the chosen path. We 
would like to create a construction that forces a choice 
between only the points in the cylinder interiors, and 
ensure all the corner points will be visited regardless of 
which path is chosen. To accomplish this, after the last 
corner of the gadget, we can have P loop back around 
along the outer edge, covering all the corner points with- 
out covering any of the interior points. Figure [2c] demon- 
strates this configuration. 

Figure l4"al shows an example usage of the gadget. The 
points in the cylinder interiors represent the clauses 
in which the variable appears. Only one set of clause 
points, either the clauses in which the positive literals 
appear or the clauses in which the negative literals ap- 
pear, can be reached. However, all the corner points 
will always be covered. The full construction will in- 
clude one of these gadgets for every variable, with each 
one passing through the points corresponding to clauses 
containing the variable's positive and negative literals. 
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Figure 3: An a-corner for a = 7r/4, with various prop- 
erties highlighted. ZABC = ir - a, ZJHL = n - a/2 

5.2 a-Corners 

We now give the full specifications of the corner con- 
structs, granting them the flexibility to bend at an arbi- 
trary angle a. Each corner consists of two components 
of P with four and three segments respectively, as well 
as four points of S. 

Figure [3] shows the full a-corner construction. The 
two components of P are A, B, C, D, E, which we refer 
to as the forward path, and F, G, H, I, which we refer 
to as the return path. The four points of S are G, H, K, 
and L. The line lengths are chosen so that the following 
properties hold: 

• \\BL\\ = \\DK\\ = e 

• AB (resp. DE) and the parallel line through K 
(resp. L) are separated by exactly e 

• G, H, K, and L are all collinear 

• The line through midpoints of A J and EE passes 
through G. 

The last property is enforced so that any points of S 
in cylinders before AB are not visible to any points in 
cylinders after DE, thus ensuring that external points 
cannot break the properties of the corner. The compact- 
ness of a-corners can be shown using simple geometry; 
an infinite strip along AB of thickness 14e is sufficient 
to contain all points of the structure for any a. 

5.3 Construction 

Using this gadget, we can now build the final construc- 
tion. 

We begin by adding an initial set of points to S which 
we refer to as "clause points" C\, . . . , G„, one for each 
clause in $. We position these points so that they lie, 
equally spaced, on a circle of radius n 2 e. This radius 
ensures that, if the lines CiCj and CkCi are parallel, 
they are separated by at a distance of least 14e. Let 



the clause strip about clauses i and j denote the set of 
points within 7s of the line CiCj. 

We now proceed to create P, adding more points to 
S as needed. As we place components of P, we require 
that all joints and a-corners be placed entirely outside 
all strips about all clause pairs, so as not to block fu- 
ture pieces. Furthermore, we require that new pieces of 
P be placed so that their cylinders do not intersect the 
convex hull of all previously placed points in S. This 
ensures that previously placed pieces do not create un- 
intended "shortcuts" that could break the properties of 
the new pieces. The exception to both these rules is the 
segment that passes directly through two clause points 
i and j, which, of course, must pierce the convex hull. 
Its adjacent a-corners will lie entirely inside the strip 
about CiCj, but must be placed outside all other clause 
strips. Note that this is always possible; beyond An 2 e 
units from the center of the clause ring, no clause strip 
intersects any other. Strips of different angles will grow 
further and further apart, creating regions of arbitrary 
size between them. 

So long as the requirements in the preceding para- 
graph are met, the start point of P can be placed ar- 
bitrarily. We then perform the following procedure for 
each variable in $, building the construction incre- 
mentally. Let x and y be the clauses in which the pos- 
itive literals of Vi occur, and z and w be the clauses in 
which the negative literals occur. We begin by position- 
ing an a-corner so that the extension of the last segment 
of the forward path passes through C x and C y . An ex- 
tra point, which we refer to as a split point, is added 
on the boundary of the first forward path segment to 
induce the splitting of the two possible paths. From 
there, both the forward and return paths are extended 
through the clause ring, with the forward path crossing 
through C x and C v . 

On the opposite side, outside the convex hull of all 
points in S so far, another a-corner is added, bending 
the path toward C Z C W . More a-corners, all bending in 
the same direction, are added as needed until one can 
be placed such that the forward path passes through G z 
and C w . Note that there must be an odd number of 
a-corners in order to ensure that C x , Cy and C w , C z 
are reachable on different curve possibilities. Once the 
paths have been extended through the clause ring and 
outside the convex hull, another split point is added on 
the boundary to collapse the curve possibilities. Finally, 
the forward path is linked to the return path, and the 
joint is added to S. At the end of the return path, more 
segments of P are added, with each joint being added to 
S, in order to move to the next variable's clause strips. 

Once this process has been completed for all variables, 
the construction is complete. Note the units in our con- 
struction are all in terms of e, so e can be chosen arbi- 
trarily. Figure I4bl shows a completed construction for 
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(a) A partial construction for formula with 12 clauses, 
showing the gadget for a single variable. The only two 
valid paths visit the either positive literal's clauses or the 
negative literal's clauses. 



(b) A completed construction for the formula $ = (iV 
yVz)A(x\/y\/z)A(xVyVz)A(x\/yVz). The upper right 
clause point represents the first clause, and the second, 
third, and fourth follow counterclockwise. 



Figure 4: Example constructions 



a very simple formula, while Figure 0a| shows a partial 
construction for a more complex formula. 

6 Result 

Lemma 3 There exists a polygonal path Q on S with 
5p(P, Q) < £ that visits every point in S if and only if 
$ is satisfiable. 

Proof. For the forward direction, assume $ has a satis- 
fying assignment. It is easy to see that our construction 
always has a polygonal path Q on S with 5f(P, Q) < £ 
that will visit every non-clause point; a-corners are con- 
structed specifically to ensure this. If $ has a satisfying 
assignment, then one of the two path possibilities in 
each variable construct will cover the the clause points 
corresponding to the clauses satisfied by that variable, 
resulting in all clause points being visited as well. 

For the backward direction, let Q be a complete polyg- 
onal path Q on S with Sp(P,Q) < e. By constructing 
each variable construct completely outside the convex 
hull of all previously placed points of S, we have ensured 
that any Q with 5f(P, Q) < £ must follow the path we 
have laid out. Each variable construct forces a choice 
between two paths, representing a true or false value 
for that variable. Since each Q visits each clause point, 



the path taken in each variable construct represents an 
assignment to the variables that satisfies $. □ 

It is straightforward to show that five a-corners is 
sufficient to move between any two strips. Thus, the 
construction is clearly of polynomial size. This, together 
with the fact that the problem is in NP, leads to the final 
result. 

Theorem 4 The Non-unique All-points Continuous 
CPSM Problem is NP-complete. 

In the construction, the only points that occur more 
than once are the clause points and the inner a-corner 
points. In all occurrences of both cases, the next point 
is always reachable from the previous point. Thus, for 
this class of problem instances, any solution to the Non- 
unique version of this problem can be converted to a 
solution to the Unique version by simply skipping the 
points that have already been visited. This shows that 
the same reduction applies to the Unique version. 

Corollary 5 The Unique All-points Continuous CPSM 
Problem is NP-complete. 

7 Restricted Problem 

The hardness of the problem stems from the fact that 
P could come within e of a point in S multiple times. 
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Thus, it is natural to ask if the hardness remains if we 
make a restriction that prevents such a situation. In 
the following section, we show that the Non-unique All- 
Points Continuous CPSM problem is polynomial-time 
solvable under the condition that, for all s £ S, the set 
{t £ [0, 1] | \\P(t), s\\ < e} is connected. 

7.1 Algorithm Outline 

Let Pi be the ith. segment of P, and let d be the cylinder 
of radius e around Pj. For simplicity, let Co = B(P(0), e) 
and Cn+i = B(P(l),e). Let Si = d H S. For a point 
s 6 S, let l(s) be the earliest occurring point of P that 
is within e of s, and let r(s) be the latest. By the 
restriction imposed in the preceding section, l(s) and 
r(s) are uniquely defined. 

A obvious preprocessing step is to confirm that all 
points of S are in some cylinder Cj. Another is to con- 
firm that S*o and SVi+i are nonempty. Since instances 
that do not satisfy these properties can be immediately 
ruled out, we assume them to be true for the remainder 
of this section. 

The ith and jth segments of P are said to be con- 
nectable via (s, t), where s £ Si and afG Sj, when the 
following properties hold. Note that s and t could be 
the same point. 

1- \J i<k< jS k -S i -Sj = ® 

2. S F (P',lt) < e, where P' is the subcurve of P from 
r(s) to l(t) 

3. VdgS.U Sj, l(v) ^ r(s) or l(t) ^ r{v) 

Let G be a directed graph whose vertices correspond 
to the segments of P, with an edge between any two 
connectable segments of increasing index. 

Theorem 6 There is a path in G from P\ to P n if and 
only if there exists a polygonal curve Q whose vertices 
are exactly S with Sf(P,Q)- 

Proof. (=>) To show the forward direction of this state- 
ment, we will construct a curve with the requisite prop- 
erties under the assumption that G has such a path, 
denoted by 1 = a\ < 02 < . . . < a m = n. Let U and s$ 
be the points in S ai connecting P a . to P ai _ 1 and P ai+1 
respectively. In other words, if C ai _ 1 and C ai are con- 
nectable via (x,y), and C ai and C Qi+1 are connectable 
via (z,w), then ti = y and s, = z. Let t\ (resp. s m ) be 
an arbitrary point in Sq (resp. S n +i); these points will 
be the first and last points of Q. 

Repeat the following for each i from 1 to to. Add ti to 
Q. Then, visit every point in the set {v £ S ai \ l(ti) ^ 
r(v) and l(v) ^ r ( s i)} ~ { s ii^}i i n order monotonic 
along the direction of P ai . Finally, add Sj to Q. Once 
this process has been performed for each i, the path is 
complete. 



We must show now that the constructed curve (i) vis- 
its every point in S and (ii) has Frechet distance at most 
e from P. 

We prove (i) first. For a given iteration i of the above 
process, assume some point in v 6 S ai where not did 
not satisfy the criteria l(ti) ^ r(v) and l(v) ^ f*(s,). 
Without loss of generality, assume v did not satisfy 
the former of the two conditions. Then by Property 
3, l(v) < r(si-i). If l(U-\) ^ r(v) as well, then v would 
have been added to Q in the i — 1 iteration. If not, then 
l(v) < r(sj_2), again by Property 3. Since l(t\) ^ r(v) 
for all v G S, it follows that v must have been added to Q 
during some previous iteration. An identical argument 
can be used to show that, if v had failed the latter con- 
dition, v would definitely be added during some future 
iteration. Thus, the final constructed curve contains ev- 
ery point in Uiefi ml wmcn by Property 1 above, 
equals S. 

Now we prove (ii). To accomplish this, we need only 
show that 6f(P',Q') < e, where P 1 the subcurve of P 
from l(ti) to r(si) and Q' is the subcurve of Q consisting 
of all points added during the ith iteration. Combining 
this result with Property 2 will prove the statement by 
concatenation. 

Fix an iteration i, and let the sequence t{ = 
vi, 1)2, ■ ■ ■ , v m — Si denote the sequence of points added 
during the zth iteration. For j £ [2, m — 1], let pj de- 
note the point on the segment P' n P ai closest to Vj. 
Let p\ — l(ti) and p m = r(si). It follows from Property 
3 and from the fact all points other than ti and Si are 
added to Q in order monotonic along the direction of 
P ai that pi ^ P2 r< . . . r< Pm- Thus, the concatenation 
of pi through p m is exactly P'. Since HujjPjH < £ for 
all j, either by definition or by virtue of being in C ai , it 
follows by concatenation that 6f(P', Q') < e. □ 

Proof. (<=) We now prove the backward direction of 
Theorem [5] Assume there exists a polygonal curve Q 
whose vertices are exactly S with 8f{P,Q) < £■ Then 
there exists a sequence of edges in Q that correspond to 
a path in G from Pi to P n . 

Let gi, . . . , q m be the sequence of vertices of Q. From 
the definition of Frechet distance, there exist continu- 
ous non-decreasing surjective functions a and r such 
that max tg [ ,i] l|P( tT (*))i Q( T (t))\\ < e - Let ti be the 
minimum value such that Q(r(f,)) = qt- As a special 
case, we impose that t m = 1. Then, let be the index 
of the segment on which P(a(ti)) lies. Note that the 
sequence A — a\, . . . ,a m is nondecreasing. We claim 
that, if ai < flj+i, then P ai and P Qi+1 are connectable 
via (q.i, qi+i). Thus, the maximal strictly increasing sub- 
sequence of A corresponds to a path in G. 

To show this, let i be such that cij < dj+i Since 
P(a(ti)) £ P ai and Q(r(ti)) — qi are within e of each 
other, it must be that £ S ai , and likewise for i — 1. We 
now show that the three properties of connectable seg- 
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ments are satisfied. The second property is immediate, 
as it is satisfied for all adjacent vertices of Q. 

For the first property, assume the set (J a <k<a S k — 
S ai — S ai+1 were nonempty, and let qj be a point in this 
set. Since qj is outside the cylinders C ai and C 04+1 , qj 
must be separated by more than e from either segment. 
Furthermore, by our restriction on the input, ^ cannot 
be in any cylinder before C ai or after C ai+1 . Thus, it 
must be that a* < Oj < Oj+i. This contradicts the fact 
that A is nondecreasing. 

For the third property, assume there was some point 
qj G Si U Sj such that r(qj) -< l{qj) and r((7 J ) -< l(qi + i). 
The former condition implies that qj must be visited 
after qi, and the latter condition implies that qj must 
be visited before qt+i- This is a contradiction. □ 

With this property in hand, a polynomial time algo- 
rithm follows by simply finding the connecting edges, 
constructing the graph G, and then constructing Q as 
in the proof above if G has a valid path. 

8 Conclusion and Open Problems 

We have shown that both the Unique and Non-unique 
versions of the All-Points Continuous CPSM problem 
are NP-complete. Furthermore, for Non-unique case, 
we have shown that a modified version of the problem 
can be solved in polynomial time. 

The Unique Subset Continuous version remains an 
open problem. However, given that the Discrete ver- 
sion is NP-complete, and that the Continuous versions 
tend to be harder, this version is almost certainly NP- 
complete as well. 
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